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A comprehensive study of the effect of wall heating or cooling on the linear, transient 
and secondary growth of instability in channel flow is conducted. The effect of viscosity 
stratification, heat diffusivity and of buoyancy are estimated separately, with some un- 
expected results. ^From linear stability results, it has been accepted that heat diffusivity 
does not affect stability. However, we show that realistic Prandtl numbers cause a tran- 
sient growth of disturbances that is an order of magnitude higher than at zero Prandtl 
number. Buoyancy, even at fairly low levels, gives rise to high levels of subcritical energy 
growth. Unusually for transient growth, both of these are spanwise-independent and not 
in the form of streamwise vortices. At moderate Grashof numbers, exponential growth 
dominates, with distinct Rayleigh-Benard and Poiseuille modes for Grashof numbers upto 
~ 25000, which merge thereafter. Wall heating has a converse effect on the secondary 
instability compared to the primary, destabilising significantly when viscosity decreases 
towards the wall. It is hoped that the work will motivate experimental and numerical 
efforts to understand the role of wall heating in the control of channel and pipe flows. 
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1. Introduction 



One of the well-known methods for delaying a transition to turbulence, for example in 
boundary layers, has been to reduce the viscosity at the wall. Such a reduction could be 
brought about by heating or cooling the surface, for example. The objective of this paper 
is to study the effect of wall heating on the instability of a channel flow. It is shown that 
heat can have surprising effects on the different mechanisms of transition. We restrict 
ourselves here to routes based on the linear eigenmodes, a direct nonlinear interaction will 
be studied in future. The emphasis here is on delaying/advancing the onset of transition 
to turbulence, rather than drag reduction in full turbulence, as achieved by adding small 
quantities of polymer. 

The critical Reynolds number for linear instability in a plane Poiseuille flow is 5772.22 



Qrszad ||l97l|)]. However, experiments usual ly find fully developed turbulence at a much 



lower R eynolds number, arou n d 1500 [see e.g . 



(1967|); 



Patel fc HeadI l|l969|) : 



Kao fc Park 



Davies Whitel ill 92Sl) : 



Naravanan fc Naravanan 



(1970|)]. It is clear that routes to turbulence 



other than the traditional Tollmien-Schlichting (TS) mechanism are in operation. The 
background noise in the flow has a major influence in delaying /hast ening transition to tur 



bulence, as well as in deciding which mechanism will be dominant 



Morkovin fc Reshotko 



1 198i<) ]. At extremely low levels of noise a traditional TS mechanism and/or secondary 



instability is likely to be followed. At intermediate levels, a tran sient growth of distur- 



bance s is the more likely mechanism for i n itial disturbance gro wth 



(2001|); 



Meseeued ()2002|) : 



Foster! (|l997|); 



Schmid fc Henningson 



Corbett fc Bottard l|2001|) ]. Once disturbance 



growth is triggered by a linear mechanism, nonlinearities are required to achieve a new 
self-sustained state. Alternatively, at higher levels of backgr o und noise, nonlinea r mech - 



anisms can directly come into play, see e.g. 



Faisst fc Eckhardd ^QQ^; 



Hof et al. 
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1 20(14^ . At present it is not un derstood exactly which route wiU be followed 



when (for a recent review on pipe flow see 



Kerswel. 



(20051)) 



The effect of wall heating on linear stability alone has been studied by several re- 
searchers. Here too, the effect of buoyancy has not been clearly quantified. To our knowl- 
edge, a detailed study of other mechanisms has not been done. Two related studies of 
transient growth had different emphasis fro m the present work. Tra nsient growth in two- 



fluid flow was studied in two-dimensions by 



understanding the effect of the interface. 



Malik fc Hoopeil (120051) with the objective of 



Biau fc Bottard l|2004(l studied transient growth 



with stable thermal stratification and concluded that such stratification is a viable strat- 
egy to control transitional flows. A more detailed retrospective on earlier work is included 
in the relevant sections later in the text. 

We consider two types of heating. The first is asymmetric, with the two walls main- 
tained at different constant temperatures. The second is symmetric, with the walls at 
one temperature and the fluid at another. Our results may be summarised as follows: 
the linear stability results are in line with the flndings of earlier studies. A decrease in 
viscosity as one approaches a wall has a large stabilising effect and vice-versa. The ef- 
fect on the line ar eigenmode of reduc ed heat diffusion (increasing Prandtl number) is 



extremely small [Wall fc WilsonI (|l996tl ] . Buoyancy has no effect up to a Grashof number 
of about 3000 and is enormously stabilising or destabilising thereafter, depending on the 
sign of the temperature difference. The Rayleigh-Benard and ToUmien-Schlichting modes 
are distinct at low Grashof number and merge at high Grashof numbers. 

The effect of heat on transient growth of instabilities is unexpected. Viscosity strat- 
ification, which is the chief player in linear instability, has no discernible effect on this 
mechanism. Increasing Prandtl number, on the other hand, has an order of magnitude 



destabilising effect. The assumption that Prandtl number may be neglected in stabil- 
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Figure 1. Schematic diagram of the chan- 
nel, asymmetric heating. The two walls are 
held at different temperatures, T^ot and 
Tcoid, SO the mean temperature profile is 
linear. 




Figure 2. Schematic diagram of the chan- 
nel, symmetric heating. Both walls are 
maintained at the same temperature, dif- 
ferent from that of the fluid. 



ity analyses is therefore conipletely incorrect for this mechanism. Secondary instabilities 
of the ToUmien-Schlichting modes are usually taken to be unimportant for channel-flow 
transition, but we find that viscosity-stratification can have a destabilising effect on these 
modes, which may make them noticeable at large temperature differences. 



2. Basic Velocity Profiles 

Two types of temperature variation, which we shall refer to as the asymmetrically 
and symmetrically heated cases respectively, are considered, see figures ^ and El These 
provide a fair sample of the type of stratification we may come across. In the first, the two 
walls of the channel are maintained at different temperatures, Thot and Tcoid- At steady 
state, the temperature within the channel varies linearly between the two. Note that for 
the unstable Poiseuille-Rayleigh-Benard configuration, the temperature difference AT 
between the bottom and top walls (and hence the corresponding Richardson and Grashof 
numbers, defined later) is taken to be positive. The sign of AT is unimportant when 
buoyancy is neglected. 

In the second case, the walls are both maintained at the same temperature, while the 
incoming fluid is at a different temperature. Such a flow is non-parallel, since the fluid 
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temperature downstream tends to equilibrate with the waU temperature, but for large 

Peclet numbers, the change in the downstream direction is very slow, and the flow may 

safely be assumed to be locally parallel and the local temperature profile to be parabolic. 

Stability analyses for more realistic temperature profiles have been conducted, without 

any qualitative difference. The temperature-dependence of the viscosity is described by 

the Arrhenius model, which works fairly well for most common liquids like water and 

alcohol. 

/i(T) = Ciexp(C2/T), (2.1) 

where Ci and C2 are constants associated with the fluid under consideration, which is 
taken in the present computations to be water. The streamwise direction is denoted as 
X, the coordinate y is normal to the wall, and z is the spanwise direction. The mean 
a:-momentum equation for a plane parallel channel flow reduces to 

ifiUy = ^i?e, (2.2) 
ax 

where the primes denote differentiation with respect to y, Re is the Reynolds number 
defined as Re = Umaxh/vref , h is the half-channel width, i^ref is the reference kinematic 
viscosity, and dP/dx is the constant pressure gradient. The viscosity ratio is defined as 
m = fJ.coid/ fJ-hot for the asymmetric heated case and m — fJ-wau/fJ-c, where c stands for 
centerline, for the symmetric case. Knowing fJ.(T) and T{y), equation 12.21 is integrated 
twice by a fourth-order Runge-Kutta method to get U. Figure 13 shows typical viscosity 
and velocity profiles for asymmetric heating. Viscosity profiles for the symmetric case are 
shown in figure 0] Corresponding velocity profiles and their second derivatives are shown 
in figure El 



Unless otherwise specified, we define the Reynolds number in terms of the average 
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Figure 3. Variation of (a) viscosity and 
(b) velocity with asymmetric heating. The 
velocity is scaled by its maximum and the 
viscosity is scaled here by its value at the 
hot wall. 
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Figure 4. Variation of viscosity, scaled 
here by its value at the centerline, for sym- 
metric heating. 
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Figure 5. Variation of non-dimensional velocity and its second derivative U" with symmetric 

heating. 



viscosity across the channel, as follows. 

Re 



Um.n.T. ^ 



(2.3) 



We can see from Table ^ that the average viscosity varies significantly with increasing 
AT, so defining the Reynolds number as above is approriate for making comparisons 
at a given Reynolds number between heated and unheated flows. On the other hand, 
it is seen from the table that the average velocity is practically unchanged by viscosity 
stratification, so the maximum velocity is a good enough velocity scale. 
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AT^K 


m 




Uavg /Uavg,m=1.0 


50 


0.39 


0.748 


1.068 


10 


0.81 


0.933 


1.006 


0.0 


1.00 


1.000 


1.000 


-10 


1.23 


1.070 


0.995 


-50 


2.51 


1.381 


0.9944 



Table 1. The dependence of the average viscosity /lavg on the viscosity ratio, m, symmetric 
heating. The last column shows the ratio of the average mean velocity to its value in the unheated 
case. The ratio is close to 1.0 in all cases. 



3. Linear stability 

3.1. The stability equations and their solution 
The disturbance quantities in normal mode form are given as 

v,fi,f = v{y),ri{y),f{y) exp[i(aa; + /3z - wi)], 



(3.1) 



where v and fj respectively are the components of disturbance velocity and vorticity in 
the direction normal to the wall, T is the disturbance temperature, a and /3 are the wave 
numbers in the streamwise and spanwise directions respectively, and w is the complex 
frequency of the wave. The linear stability equations may be derived to be 

ia [{v" - {a^ + /?2))(J7 - c) - U"v\ = ^[l^ [v'" - 2{a^ + P>" + + H^v] 
+ ^T'2[v"' - ia' + P')v'] + ^T"[v" + (a' + p')v] + ^T"[v" + {a' + p')v] 



+ -r^U'T'T——iaT 
dT^ Re 



dT-^ 



dT^' 

+ ^[U'f" + 2U"f' + (a'U' + U"')f] + 2^U'T'f' + ^T"U'f 



(3.2) 



a{c - U)r, + ipU'v = [^,[r," - {a' + (3')r,] + ^T'rj' - il3^{U"f + U'f) 



dT' 
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{U - c)f - T'v 



(3.3) 



(3.4) 



RePr 

where c = uj/a. Equations 13 . 21 and 13 .31 respectively are the Orr-Sommerfeld and Squires 
equations, modified here to account for the effects of viscosity variations, temperature 
fluctuations and of buoyancy. A Boussinesq approximation has been made. If buoyancy 
were to be neglected, the above equations would be equivalent to those of Wall & Wilson 
(1996). 

The Prandtl number is defined as Pr = v/k where k, is the coefficient of thermal 
diffusivity. The Grashof number is Gr = g^ATh^ /v"^, 7 being the volume coefficient of 
expansion. Eauations l3.2l to l3.4l form an eigenvalue problem with the boundary conditions 



?;(±1) = v'{±l) ^ 77(±1) = f (±1) = 0, 



(3.5) 



and are solved using a Chebyshev collocation spectral method. We perform a temporal 
stability analysis, the growth rate of the disturbance is obtained from the imaginary part 
of c. 



3.2. Effect of viscosity variation 

In order to isolate the effect of viscosity variation, the Prandtl number and the Grashof 
number are set to zero. With symmetric heating, we expect that if the viscosity decreases 
as one approaches the wall, the fuller velocity profile will result in a stabler flow. This 
expectation is realised, as seen from the neutral stability boundaries in figure 

However, in a channel fiow where one wall is maintained at a constant high temperature 
and the other wall is kept cold, the viscosity decreases towards one wall and increases 
towards the other. It is not a priori e v ident what the effect on the linear stability will 



be. It was found by (Potter fc Graben l|1972|) that any temperature difference between 
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the walls is always destabilising. However, 



IWa,]] WilsonI 



(122a) found, using four dif- 



ferent viscosity models, that a temperature difference almost always stabilises the flow. 
The apparent contradiction is because the former work compared results for heated and 
unheated flow maintaining the input power constant, while the latter made comparisons 
at a given Reynolds number. Since the flow rate for a given input power is higher for the 
heated case, the resulting Reynolds number is higher. The sta b ility o f viscosity-stratified 



chann el flows was also studied by 



Pinarbasi fc LiakoDoulof 



19931) with siniilar conclus i ons and of boundary layer flows by 



Wazzan et al. 




196 




197' 



Strazisar et al 



(1977|); 



Cra 



Schafer et al 



Schafer 



iki lll97ll) : 




Kao 



Herwie 




1963); 



199 



In the present paper, we define the Reynolds number in terms of a verage viscosity, 



and c ompare results at a given Reynolds number. In agreement with 



Wall fc Wilson 



1 1996^ . for asymmetric heating, we find that any temperature difference is stabilising, in 



terms of the le ast stable (two-dimensional) linear mode (figure (Tj) . We have confirmed 



Sameen 



1 20041) ] that the production of disturbance kinetic energy is reduced at the 
cold wall and increased at the hot wall compared to the unheated case. The dissipation 
is similar in all cases. The highly oblique modes, unlike the two-dimensional ones, are 
practically unaffected (not shown). This observation will assume significance when we 
discuss transient growth. 

3.3. Effect of heat diffusivity 

We know that for liquids such as water, heat diffuses slower than momentum, so the 
assumption of Pr = is not justifiable. Surprisingly however, the linear stability, as 
measured by the least stable eige nmode, is practically unaffected by a decrease in heat 



diffusivity 



Wall fc WilsonI l|1996|) ]. Present computations confirm this (figure (SI. How- 



ever, the prevalent conclusion that heat diffusivity does not affect fiow stability, and 
therefore that the Peclet number may be set to zero in stability analyses, is shown in the 
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Figure 6. Effect of viscosity variation on linear stability, symmetric heating. 
For unstratified flow, i.e., at m = 1.0, Rccr is 5772.2. 
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Figure 7. Stability boundaries for various viscosity ratios. Asymmetric heating, 

Tcoid = 295°K. 

next section to be incorrect. Increasing the Prandtl number to 0(1) vaules can enhance 
transient growth by an order of magnitude. 

3.4. Effect of buoyancy: the Poiseuille-Rayleigh-Benard problem 

We consider the asymmetrically heated case here. When the upper wall is cold relative 
to the lower one, the resulting unstable stratification of density leads, at low flow rates. 
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Figure 8. Most unstable eigenvalue at various Prandtl numbers for different 
AT at Q = 0.9 and Re = 1000, 5000. The effect of Prandtl number is negligible. 



to a buqvancy driven instability s imilar to the Rayleigh-Benard 



Turned l)l959|) : 



been studied W 



is reviewed in 



atten fc Leg 



Zhang_^i_a^ 



Flatten fe Legro! 



998) 



(1£SJ) 



Chandrasekhar 




ros _ 1984 )] . The eff e ct of mean shear on this instability has 



Deardorfi J 196E K for ex ample, and this problem 



Mahaian et al 



1 198a) . On the other hand, the 



effect of buoy ancy on the ToUmien - Schlichting; modes in plane Poi seuill e flow has been in- 



vestig ated by 



Cage ReidNl 9681) : 



TveitereidI l|l974l) and 



Fuiimura fc Kelly 



1 19881) . Several approximations were made in these early studies. For example, viscos- 
ity variations were neglected and the base flow was taken to be parabolic. It was found 
that a critical Reynolds number always exists for any level of density stabilisation, while 
there exists a Richardson number {Ri = Gr/Re^) above which the flow is stable for all 
Reynolds number. 

In figure El the critical Reynolds number, i?ecr for a temperature differ ence of 25°K 



is plot ted f or various Richard son numbers. The trends are the same as in 



{mm and 



Gage fc Reic 



TveitereidI 1)19741 ). but there are minor numerical discrepancies, which we 



attribute to the more appropriate velocity and viscosity profiles used here. The effect of 
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Figure 9. Neutral stability Reynolds number as a function of Richardson number, 
Pr = 1.0. Note that for stable stratification (AT < 0), the Richardson number is 
negative, its absolute value is plotted here. The portion to the right of the neutral 
points is unstable. 

buoyancy are negligible when the Richardson number is below 10"*, and of either sign. 
At higher Richardson numbers, for unstable stratification, figure 1^1 shows that the flow 
is highly destabilised. The stability boundaries are plotted in figure ^1 in terms of the 
Grashof number, a given Grashof number being more simple to achieve experimentally. 
Distinct modes of Rayleigh-Benard type and of ToUmien-Schlichting type are evident 
at intermediate levels of Gr. The modes merge at Grashof numbers above ^ 25000. 
The numerical value at the bifurcation point varies slightly with Prandtl number and 
temperature difference. The unstable region in the Grashof-Reynolds parameter space is 
shown in figure ITTI 

We have discussed the stability in terms of the most unstable linear mode. However, a 
transient growth of decaying modes can often be the dominant mechanism of transition 
to turbulence in channel fiows. We shall see in the next section that the effect of heat on 
flow instability throws up several surprises. 
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10000 



Figure 10. The neutral curves for unstably thermal stratified flow at Prandtl 
number 7.0, AT = 25°K. The second (Rayleigh-Benard-like) mode starts appear- 
ing at Gr = 2695 and merges with the Tollmien-Schlichting mode at Gr = 25186. 
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Figure 11. A consolidated picture of the variation of the critical Reynolds number 
with the Grashof number. The Tollmien-Schlichting mode is shown by the solid 
line, and the Rayleigh-Benard mode by the dashed line. The region above the 
curves is unstable. 

4. Transient Growth 

The linear stability operator is not self- adjoint, and the resulting non-orthogonality of 
the eigenfunctions is known to be able to give rise to large levels of transient growth of dis- 
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turbance kinetic energy even when all individual eigenmodes are stable. In wall-bounded 
flows, transient gro wth is mainly caused by the i nteraction between the O rr-Sommerfeld 



and Squire modes |Reddv fc HenningsonI l|1993(l 



Criminale et al 



(2003)] from the cou- 



pling term, —i/SU', appearing in Squire's equ ation. The mo s t likel y structures arising due 



to transient growth are streamwise streak s 




Reddv et al 



(1993); 



Reddv fc Henningson 



Schmid fc HenningsonI ()200ll) ]. We use the standard approach for comput- 



ing the maximum transient growth. 

The effect of viscosi t y-stra tification, in contexts other than heat 



(20051); 



Chikkadi et al 



Malik fc Hoopei 



( 2005|) ] has been addressed before, though not completely. The 



effect of buoyancy has been studied under stable stratification alone by 



Biau fc Bottaro 



1 200^. The effect of heat diffusivity on transient growth has not been studied before, to 
our knowledge. 

The disturbance kinetic energy, g{t) 



9{t) 



Schmid fc HenningsonI (|200l|) ]. is written as 



II 



\e 



Mmi \\<mi 



(4.1) 



Its time evolution is represented by the matrix dn/ dt = — zAk, where n — i^n)'^ 
and A = diag{uji,uj2, ---^n}, where Kj is the j*'' expansion coefhcient of the eigenfunc- 
tions of the linear modes which are the dominant contributors here. The superscript T 
denotes transpose. Maximising equation (|4.1|) for all possible initial conditions k(0), 



G{t) = max g(i). 



(4.2) 



We then define Gmax as the maximum over time of G{t) for one particular Re, a and 
AT, see figure |ni 

The contour plot for G,nax for unheated flow is shown at Re — 1000 in figure IT^ for 
comparison with the results for heated flo w to follow. A maximum grow th of Gmax = 196 



is obtained for a = 0.0 and /3 = 2.05 [see 



Reddv fc HenningsonI (|1993(1 ] 
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Figure 12. The energy amplification evolution for various viscosity ratios for 
Re — 3000, Q = 1, asymmetric heating. The thick dashed curve is for unstratified 
flow, and the solid lines are for m = 1.2, 1.65 and 2.5 in the order indicated in the 
figure. 
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Figure 13. The contour of Gmax (the maximum over time of G(t)) for Re = 1000 



in the a — P plane, AT = Q° K. This matches well with 
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m 

Figure 14. The variation of Gmax at a = and P = 2 for various viscosity ratios 
at Re = 1000, asymmetric heating. The maximum deviation of Gmax from the 
unheated value of 196 is only 3%. 

4.1. Effect of viscosity stratification 

As before, we first take the Prandtl number to be zero, i.e., assume that temperature 
fluctuations diffuse away instantaneously. We also neglect buoyancy, in order to isolate 
the effect of viscosity stratification alone. For the asymmetrically heated case, the growth 
of kinetic energy is seen in figure lT^ to change only marginally with heating. The example 
shown in figure IT^ is for Re = 3000 and a — 1, but viscosity stratification has very little 
effect on transient growth at any value of Re and a. The effect of asymmetric heating 
is quantified in figure IT^ in terms of Gmax at a = and (3 = 2. There is a marginal 
stabilisation with viscosity stratification. This result is in line with the result for linear 
stability, but much smaller in magnitude. The contours of Gmax for symmetric heating 
are plotted in figures 1151 in the a ~ (3 plane, and the variation of Gmax with viscosity 
ratio is plotted in figure [TCI at a = 0. There is again a slight stabilisation with increase 
in viscosity stratification. 

The insignificant effect of viscosity stratification is consistent with our recent study 
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Figure 15. The contour of Gmax for symmetric heating at Re = 1000 again in 
the a — P plane, m = 0.4. 
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Figure 16. The Gmax variation at a = 0.0 and 13 = 2.0 for various viscosity 
ratios at Re = 1000, symmetric heating. The maximum deviation of Gmax from 
that for unstratified flow is only 13%. 



Chikkadi et al 



(2003)]. As 



of transient growth in two-fluid and non-Newtonian flows 
discussed there, the U" term, which affects the least stable eigenmode dramatically, has 
no effect on streamwise vortices arising from a = 0, which dictate transient growth. The 
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Figure 17. (a) The eigenspectra for two extremes of viscosity stratification, 
m — 0.4 and 2.5, for a = 0, i?e = 1000 and /3 = 2.0, symmetric heating, (b) 
The corresponding eigenfunctions of the first two unstable eigenvalues in each 



eigenspectrum, and typical eigenfunctions at a = are shown in figure 1171 to be very 
similar at two extremes of viscosity stratification. Eauation l3.3l drivcs the dynamics rather 
than equation 13.21 under these conditions, and the terms containing viscosity gradients 
have been verified numerically to be small. 

4.2. Effect of heat diffusivity 

It has been seen that the Prandtl number has a marginal effect on the most unstable linear 
mode. In contrast, we find here that reducing heat diffusivity has a large destabilising 
effect on the transient growth of disturbance kinetic energy. A dimcnsionless quantity for 
measuring growth is the energy norm defined as 



E= I \up\^ + \vp\'' + \wp\^ 



\T\^dy. 



(4.3) 



There is some flexib ility in defining the measure , but the results are no t expected to 



change qualitatively 



Hanifi fc HenningsorJ l|l998l) : 



Riau RottaTfJ 1)20041) ]. In figure EHl 



for a temperature difference of 25° K at a Reynolds number of 1000 the effect of Prandtl 
number is shown. As the Prandtl number is increased from 10^* to 1, the transient 
growth is seen to increase dramatically. The large destabilisation comes from a new two- 
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Figure 18. Asymmetric heating: contour plot of Gmax for T = 25° A', Re — 1000 
for various Prandtl numbers. Note that for Pr = 1 the scale employed is different. 

dimensional transient growth. This is true for symmetric heating as well (not shown). 
We now have a situation where transient growth dominates, but not via the standard 
streamwise streaks and streamwise vortices. 

4.3. Effect of unstable density stratification 
In their studies of stable thermal stratification 



Biau & Bottaro 




have found that 



as stratification increases flow becomes increasingly stable, both in terms of exponential 
growth as well as transient growth. Viscosity variations were not accounted for in their 
calculation. In this paper, we concentrate on unstable thermal stratification. 

Figures lT^ to l^ show contour plots of maximum growth of transient energy for various 
Grashof numbers at increasing Prandtl number for a temperature difference of AT = 
25°K. As expected, when heat diffusivity is high, buoyancy has little effect. For Pr ^ 
0.01, the Prandtl number dictates the instability, and buoyancy has very little effect. At 
Pr — 1, however, the situation is completely different. For Grashof numbers of 1000 and 
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Figure 19. The transient growth at various Grashof number for AT = 25° AT, 
Re — 1000 and Pr = 10~*. Since temperature perturbations diffuse away rapidly, 
buoyancy does not have much effect. 




a a 

Figure 20. Same as figure EH but for Pr = 10"^. The Grashof number has no 
effect upto a value of ~ 10''. At Gr = IC* a new growth appears at /3 = 0, which 
win dominate at higher Prandtl number. 

below we see extremely large levels of subcritical transient growth. This growth is two- 
dimensional. Above this Grashof number of course, a linearly unstable mode exists. The 
transient growth in unhealed channel flow is well known to display itself as streamwise- 
independent structures, like streaks and vortices. Our results indicate that such structures 
will not be much in evidence in heated flows at realistic Prandtl or Grashof numbers. 
Rather, a spanwise-independent growth occurs. Experimental and numerical veriflcation 
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Figure 21. Same as figure \W\ but for Pr — 10~^. The new subcritical mode is 
now dominant at low Grashof number as well. The spanwise-independence of the 
largest transient growth, making the Poiseuille-Rayleigh-Benard problem essen- 
tially two-dimensional, is unusual. 

of this kind of transient growth could have interesting consequences for wall heating as 
a control option. 



5. Secondary Instability 

A flow containing linear modes (either growing or decaying) of sufficient amplitude 
Ap can become unstable to new secondary modes of instability. In unhealed channel 
flow, secondary instability is considered unimportant, since it may only play a role when 
external disturbance levels are extremely small. We show here that viscosity variations 
can significantly destabilise the secondary mode, thus making it more relevant to the 



The approach here is as in 
are decomposed in the form 
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1 19881) . AU flow variables 



u{x,y,z,t) = U{y) + ApUp(x,y,t) + AsUsix.y, z,t), 



(5.1) 



where Up is the linear instability mode of the previous section. A subscript p has been 
introduced in this section to denote primary mode of instability. As is the amplitude of 
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Figure 22. Same as figure [T^ except that Pr — 1. The transient growth is even 
higher than before. The region inside the thick curve for Gr = 10^ is hnearly 
unstable. 

the secondary. Note that since only the least stable linear modes are relevant here, it is 
sufficient to consider two-dimensional primary modes, by Squire's theorem. 
The secondary perturbation quantities are assumed to be of the form 



+ (w_,w_,u._)(y,i)e'("-^-''=^) +c.> 



(5.2) 



where a+ and a_ are the wave numbers of the secondary waves in the streamwise di- 
rection, Ps is the wave number in the spanwise direction. The direct interaction between 
primary waves is assumed to be negligible. For the flow under consideration, the growth/ 
decay rates are so small that dAp/dt can be neglected during one period of time, and the 
primary flow may be taken to be periodic. Substituting these decompositions in the mo- 
mentum equations, eliminating pressure and neglecting non-linear terms in the secondary 
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disturbance, we arrive at the secondary disturbance equations. On averaging these over 



X, z and t, only the resonant modes survive, which are given by 



a-i- + Of- 



(5.3) 



The cases of a+ = a/2 and a+ = a are called the subharmonic and the fundamental 
modes respectively. Using continuity the streamwise component of secondary disturbance 
velocity is eliminated and we get the secondary perturbation equations: 



-D^ + = ~sAf+ + {AD - ia+{DU))v, 
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D — iuz 



i{a + a^)upD — ia-{Dup) + Vp(^ ^^'^ + D^^ 



fl 



(5.5) 



where A = [ia^U + — ^(P — ^' D], 



D = d/dy. Equations H5.4|l & (|5.5|) and 



complementary equations for and fl are solved using a Chebyshev collocation spectral 
method, with the boundary co nditions , Vfi , W a — at y — ±1. The dispersion relation is 
F{Ap, (3s,m, Re,a,c,) = [see^^ 



Herbert! 1)19831) ] . The growth rate is highly sensitive to the 
primary amplitude level Ap, and increases with increasin g An. The prese nt computations 



are validated b y comparing with the unstratified case in 
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1983 ) as discussed in 



5.1. Asymmetric heating 

A value of Ap — 0.01 is taken, to be representative of an intermediate level of primary 
disturbance. The variation of secondary growth rate with the spanwise wave number 
for various viscosity ratios is plotted in figures l^ As the temperature difference increases, 
a second highly unstable mode appears. This mode is three-dimensional, the spanwise 
and streamwise wavelengths being very close to each other. The modes which are closer 
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Figure 23. Asymmetric heating: dependence of growth rate on spanwise wave 
number of the secondary disturbance for various viscosity ratios, subharmonic 
case, a = 1.0, Ap = 0.01, Re = 5000. 

to two-dimensional arc now stabilised. A nonlinear or transient growth triggered by this 
new mode could mean that transition to turbulence proceeds somewhat differently, but 
further studies are needed to evaluate this. 

While we have taken the amplitude of the primary mode to be constant, it is in 
fact, at these Reynolds numbers, a known slowly decaying function of time. Integrating 
instantaneous results over long times, we can compute the time dependence of amplitude 
of the secondary mode. This approach is a counterpart in time of the assumption of 
parallel flow in flows which vary slowly in x. The amplitude of the subharmonic secondary 
mode As is shown as a function of time in figure 1^ At low initial Ap, secondary modes 
are always stable while for higher Ap significant growth is displayed up to large times. At 
low Reynolds number, the initial Ap needed for a sustained secondary instability growth 
is very high. 
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Figure 24. Asymmetric heating: variation with time of the amplitude of sec- 
ondary disturbance for two sets of initial Ap. a — 1.0, Ps = 1-0, Re = 5000, 
Subharmonic mode. 



In figurel^the phase shift Pg 



- LOg is shown as a function of the spanwise wave 



number. The phase locking of the subharmonic wave (i.e., where Pg is zero) is achieved 
at an earher (3s than for the unstratified case. 

5.2. Symmetrically heated channel 

Figure 1251 shows the secondary perturbation growth rate variation with spanwise wave 
number (3^. A stabihsation with increase in viscosity ratio, especially when m > 1, is 
evident. This behaviour is remarkable, being counter intuitive and opposite to the be- 
haviour of the primary instability mode. The phase locking behaviour is similar to the 
asymmetrically heated case. 

We have seen that in linear disturbance growth, the mean velocity p rofile (via the U" 



term) has a dominant role. In the case of secondary growth as well, 



Orszag fc Patera 



1 1983|) have argued that inviscid effects are dominant, and act through vortex stretching 
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Figure 25. Asymmetric heating: variation of Ps with spanwise wave number for 
various viscosity ratios, subharmonic case, a = 1.0, Ap = 0.01, Re = 5000. 
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Figure 26. The dependence of growth rate on spanwise wave number of the 
secondary disturbance for various viscosity ratios, subharmonic case, a = 1.0, 
Ap = 0.01, Re = 5000, symmetric heating. 
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Figure 27. The new mode of instability arises from changes in the velocity profile. 
This is verified by switching off the viscosity gradient terms in the secondary 
stability equation. The dashed lines are for the full equation and the filled circles 
are with the viscosity gradient set to zero in the instability calculations, but 
retained in the mean flow computations. AT = 50(m = 2.5), q = f .0, Qs = 0.5, 
Ap — 0.01, Re = 5000, asymmetric heating. 

and tilting. We are not able to make a conclusive statement on this, but it seems that 
heating affects secondary instabihty by an inviscid mechanism, through changes in the 
velocity profile. This is demonstrated in figure 1271 where it is seen that switching off 
the viscosity gradient terms makes little difference to the result. Also, it is not evident 
why the sign of instability is opposite to that of primary modal growth. The relevant 
inference here is that viscosity stratification alone can have unexpected effects on the 
various mechanisms leading to transition to turbulence. 
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6. Conclusion 

The control of the flow using heating or cooHng of the surface has long been prac- 
tised, especially in open flows. In plane channel flow we have conducted a comprehensive 
study of the effect of heat, and show that there is no unique direction (either towards 
or away from stabilisation) in which the flow responds. Linear stability results are con- 
sistent with earlier studies in that the most unstable linear mode is suppressed when 
viscosity decreases towards the wall. Also the effect of Prandtl number is negligible. For 
Grashof number between about 3000 and 25000, separate modes of Rayleigh-Benard and 
ToUmien-Schlichting instability are evident, the former is at low Reynolds number. At 
higher levels of buoyancy, the modes merge. 

The transient growth of disturbances is unaffected by viscosity stratification, but 
hugely increased by reduced heat diffusivity. Both of these are counter to the effect on 
the least stable linear mode. The Prandtl number is thus not an unimportant parameter, 
as was hitherto assumed. Transient growth is also very high in the presence of buoyancy 
of the appropriate sign. With increasing Prandtl and/or Grashof number, the growth is 
two-dimensional, not in streamwise streaks, which is quite unusual for transient growth. 

Secondary instability of Tollmien-Schlichting waves is not considered an important 
player in the transition to turbulence in unheated channel flow, unless the free stream is 
unrealistically quiet. We show that a new destabilising mode appears with heating, for 
the case where viscosity is decreasing towards the wall! 

It is hoped that this work with give impetus to experimental and computational studies 
to check these predictions and to explore wall heating in all its aspects as a control 
strategy for channel and pipe flows. 

We are indebted to Prof. O.N. Ramesh for discussions throughout the course of this 
work and for reading the manuscript. A part of the work was carried out during the 
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